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Abstract 

Similar to the treatment of dense gases, fluid-dynamic equations for the 
dynamics of congested vehicular traffic are derived from Enskog-like kinetic 
equations. These contain additional terms due to the anisotropic vehicle in- 
teractions. The calculations are carried out up to Navier-Stokes order. A 
linear instability analysis indicates an additional kind of instability compared 
to previous macroscopic traffic models. The relevance for describing granular 
flows is outlined. 
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An efficient infrastructure is an essential precondition for every industrialized country. 
Therefore, the considerable deterioration of the traffic situation on 'freeways' during the 
last decade is a serious problem. Not only does impeded traffic cause economic losses of 
many billions of dollars each year, it also produces serious ecological damages. Thus, great 
efforts have been made to develop methods for traffic optimization and forecasts, for which 
reliable traffic simulations are necessary. Consequently, research on traffic dynamics has 
recently become a very important topic. During the last years, numerous results have 
been published on microscopic models |l]-§J, including cellular automata models and 
molecular-dynamics-like models 0,§J, kinetic models fl-|T0|, and macroscopic (fluid-dynamic- 



like) models JTT],||-0, aiming at an understanding of stop and go traffic. The topic is related 
to the fields of non- linear dynamics ||TTJ, phase transitions [|iT^,|TT|j, and stochastic processes 



EH- 

Macroscopic traffic models are not only suitable for on-line simulations of traffic networks, 
but also for analytical investigations. Some very recent publications proposed to derive re- 
alistic traffic equations from the 'microscopic' dynamics of driver-vehicle units via a kinetic 
approach PHIT|. However, for the following reasons, these attempts have been only partly 
successful. Either the models treat the vehicles like point- like objects [§-0]. Then, the result- 
ing macroscopic models are only valid for free traffic flow. Or the models take into account 
the finite space requirements of vehicles PJ8|-[I0|. However, the calculations were carried 
out up to Euler order only, based on the (zeroth order) approximation of local equilibrium. 
This assumes that the form of the equilibrium velocity distribution remains unchanged in 
dynamic situations, but it is given by the local values of the density p(r, t) , average velocity 
V(r,t), etc. It will be demonstrated that corresponding traffic models are not even valid 
in linear approximation, since their linear stability analysis gives totally misleading results. 
To solve this problem, we must calculate the (first order) Navier-Stokes corrections of the 
macroscopic traffic equations, which take into account the structural change of the velocity 
distribution in inhomogeneous traffic situations. 

For ordinary gases, the Navier-Stokes terms (transport terms) are calculated from the 



kinetic equation by means of the Chapman-Enskog method [12|]. An approximate, but more 



intuitive method bases on the relaxation time approximation [|T3| , |i4f (see below). For kinetic 
traffic models the situation is more involved, since the interaction term does not vanish in 
equilibrium situations. Now, this problem has been solved. In the following, we will show 
how to derive realistic traffic equations which include corrections due to vehicular space re- 
quirements and Navier-Stokes terms. Since these non-trivial corrections change the structure 
of the equations, they cause an additional instability. For the same reason, the low-density 
regime does not allow an extrapolation to situations at high densities. The presented method 
is also relevant for understanding instabilities in granular flows, since granular collisions are 
also not energy conserving. Recently, a lot of publications tried to tackle these interesting 
problems with fluid-dynamic equations derived from the associated kinetic equation |15|-[T7|. 



However, most of them are restricted to Euler order or the low-density regime |T6|,|T7|, thereby 
neglecting relevant sources of stability and instability. Probably for this reason, these ap- 
proaches have not been fully successful in describing the formation of density waves in sand 
which is falling through a vertical pipe (T7|,[TB . 



The model: Since the number of vehicles on a (for simplicity: circular) freeway is con- 
served, the kinetic traffic equation for the phase space density p(r,v,t) of vehicles with ve- 
locity v = d t r at place r and time t has the form of a continuity equation with a sink/source 
term: 

d t p + d r (pv) + d v (p d t v) = (d t p) ss . (1) 

As usual, we will assume the acceleration law d t v = (Vo — v)/t, where r denotes a density- 
dependent acceleration time and Vo the desired velocity, which is assumed to be the same 
for all vehicles, here (case of a speed limit). The sink/source term (d t p) ss originates from 
sudden (non-differentiable) velocity changes. It splits up into a velocity-diffusion term due 
to fluctuations of the acceleration behavior ('imperfect driving') and an interaction term: 

(d t p) ss = d 2 v (pD) + (d t p) hlt . (2) 
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The interaction term reflects sudden deceleration processes. In analogy to the Enskog theory 
of dense gases |H]J2(J and granular media |TJj , but with an interaction law typical for vehicles 



H|[], it is of the form 

(d t pU = (l-p)x(r + l,t)B(v) (3) 
with the Boltzmann-like interaction function 

B(v) — / dw (w — v) p(r,w,t)p(r + s,v,t) (4a) 



w>v 



dw {y — w)p(r,v,t)p(r + s,w,t) . (4b) 



According to this, the phase-space density p(r, v, t) increases due to deceleration of vehicles 
with velocities w > v , which cannot overtake vehicles with velocity v. The density- dependent 
probability of immediate overtaking is represented by p. A decrease of the phase space 
density p(r, v, t) is caused by interactions of vehicles with velocity v with slower vehicles 
driving with velocities w < v. The corresponding interaction rates are proportional to 
the relative velocity \v — w\ and to the phase space densities of both interacting vehicles. 
By s(V) = Iq + l(V) (~ vehicle length + safe distance) it is taken into account that the 
distance of interacting vehicles is given by their velocity-dependent space requirements. 
These cause an increase of the interaction rate, which is described by the pair correlation 
function x( r ) = [1 — p{ r ^) s ] 1 & t the 'interaction point' r + 1. A more detailled discussion 
of the above kinetic traffic model is presented elsewhere ||J§. By describing the individual 
acceleration behavior via dtv = (Vq —v)/t and by introducing a velocity-diffusion term as 
source of velocity variations, it improves the original approach by Prigogine and Herman 
H||, which assumes a relaxation of the actual phase-space density to a desired one 0. 

Now, we will focus on the the macroscopic equations for the spatial density p(r, t) = 
J dv p(r,v,t), the average velocity V(r,t) = J dvvp(r,v,t)/p(r,t), and the velocity variance 
0(r, t) = J dv[v — V(r, t)] 2 p(r, v, t)/p(r, t). These are obtained by multiplying the kinetic 
equation with v k , integrating with respect to v, and a number of straight-forward calulations 
@f|. In order to underline the crucial results of this paper, we will first discuss the case of 



negligible space requirements (s, /, Iq <C l/p(r, t)), in which the macroscopic traffic equations 
read 

d t p + Vd r p=-pd r V, (5) 
d t V + Vd r V = -1/ pd r V + (V - V)/r - (1 - p)V , (6) 

d t e + vd r e = -2v/pd r v - i/pd r j + 2(d - e/r) 

~(1-P)J- (7) 

Here, P = p© denotes the 'pressure' and J{r,t) = p(r,t)T(r,t) = J dv[v — V(r,t)] 3 p(r,v,t) 
the flow of velocity variance. @ is the expected continuity equation for the density. In com- 
parison with the conventional Euler equations for ordinary gases, the velocity equation @ 
and the variance equation (|7|) contain two additional terms, each of which breaks momen- 
tum and energy conservation. The respective last terms result from the anisotropic vehicle 
interactions, while the previous terms reflect acceleration behavior and velocity fluctuations. 

It can be shown that the kinetic traffic equation has the Gaussian equilibrium solution 
Po(v) = p(27T@) -1 / 2 exp[— [y — V) 2 / (20)], which additionally fulfills the implicit equilibrium 
relations V — V — r(l — p)px® and 6 = Dr. If the local values p(r,t), V(r, t), and 
0(r, t) are inserted, instead, we obtain the Euler approximation. It leads to J = pT w 
Pj2l[|. However, in inhomogeneous traffic situations the form of the velocity distribution 
P(v;r,t) = p(r,v,t)/ p(r,t) changes due to the finite adaptation time r which is needed 
to reach local equilibrium. In relaxation-time approximation []5| |l3|,|l4| we find the Navier- 
Stokes correction T = — 3-\/7rO/[(l — p)p]d r Q. The corresponding instability diagram is 
depicted in Fig. 1 for the following model functions approximating empirical results j||§: 
T (p) = 8s/[0.97exp(-p/16km" 1 ) + 0.03], p{p) = exp(-p/16 km" 1 ), and D = 0.03V 2 /r(p). 

The instability diagram is obtained by (i) assuming a small periodic perturbation 
5g(r,t) = g exp[ikr + (A + iu)t] of the macroscopic traffic quantities g G {p, V, 0} rela- 
tive to the stationary and spatially homogeneous equilibrium solution g e (p) (go being the 
amplitude, k the wave number, A the growth rate, and uj the frequency of the perturbation), 
(ii) inserting g(r,t) = g c + 5g(r,t) into the macroscopic traffic equations, (iii) neglecting 
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quadratic terms in the small perturbations 8g(r, t) <C g e , (iv) determining the three complex 
eigenvalues A = A + iu of the linearized equations in dependence of p and k. Equilibrium 
traffic flow is unstable, giving rise to the well-known phenomenon of stop and go traffic p|,|TT| , 
if at least one of the growth rates is positive, i.e. max A > 0. Therefore, the instability 
diagram shows max X(k,p) if this is greater than zero, otherwise 0. 

It is interesting to compare Fig. 1 with the instability diagram in Fig. 2 which corresponds 
to the Euler approximation with T = 0. This shows that the curious maxima in the middle of 
Fig. 1 are an effect of the deformation of the Gaussian velocity distribution in inhomogeneous 
situations. Since they originate from the terms containing T, they are directly connected 
with the dynamic variance equations. Therefore, they could not be discovered in previous 
traffic models which eliminated @(r, f) by means of approximations of the kind G(r, t) pa 
e c (p(r,t),V(r,t)) H. 

Apparently, the Navier-Stokes corrections do not cause a stability of equilibrium traffic 
flow with respect to perturbations of large wave numbers k (i.e. small wave lengths i = 2ir/k). 
This surprising result is a serious problem for a numerical solution of the above equations. 
It comes from the fact that (shear) viscosity terms are missing due to the spatial one- 
dimensionality of traffic flow. This problem vanishes when corrections due to vehicular 
space requirements are taken into account (/o = l/pmax, I = 0.8 s- V). One would not expect 
this, since, for ordinary gases, the structure of the fluid-dynamic equations does not change 
at high densities. The only thing what changes are the constitutive relations for V and J 
T9|fl. However, in the case of traffic dynamics the situation is completely different. Since 



the anisotropic vehicle interactions do not fulfil momentum and energy conservation, they 
lead to contributions that cannot be absorbed by modified functions V and J . Whereas 
the vehicle density still obeys the continuity equation (J^), the structure of the velocity and 
variance equations changes considerably (7 = (1 — p)x)' 

8 t V + Vd r V = -[l/p + 7 s(l + p X l)]Q d rP 

+ 1P (2sJe/^-pQ X l 2 /V)d r V 
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- [i + 7 ps/2] a r e 

- 1 s(s/2 + pxl 2 /2)Qd 2 r p 
+ r]/pd 2 V - 1P s 2 /4d 2 Q 

+ (V -V)/r- 1P e (8) 

d t e + vd r e = -[2 + 7 ps]e <9 r v 



+ 2-fpsJQ/TTd r Q - jps 2 e/2d;V 



+ ips 2 yjQ/Tr d 2 Q + 2{D - 0/r) . (9) 

This result is valid up to Euler order. It has been obtained by evaluating the kinetic equation 
on the assumption of a Gaussian velocity distribution P (v;r,t) and second order Taylor 
expansion of the functions p and x with respect to s and I around r, thereby neglecting 
products of partial derivatives. Fig. 3 shows that the finite space requirements of vehicles 
cause the desired stability of equilibrium traffic flow with respect to perturbations of small 
wave lengths. This comes from the finite viscosity coefficient r\ = •yp 2 [s 2 \JQ/ti— p@x/ 3 / (2V)]. 

Nevertheless, the result is not even correct in linear approximation, since inhomogeneous 
traffic again changes the form of the velocity distribution. To calculate the Navier-Stokes 
corrections, we must derive an equation for the small deviation 5p(r,v,t) from po(r,v,t) = 
p(r, t)P (v] r, t) which is caused by inhomogeneities d r p, d r V, and d r Q. This has been done by 
means of the relaxation time approximation P, p!3| , ^4| which assumes that (i) the deformation 



of the local equilibrium distribution Po(v; r, t) is caused by the interaction term, (ii) the non- 
equilibrium corrections of the latter can be adiabatically approximated by — Sp(r, v , £)/t , 
where — 1/ro denotes the slowest eigenvalue of the linearized interaction operator, (iii) l/r 
is of the order of the vehicular interaction rate 

i 

1 -p)x{r + l,t)/p(r,t) 



To(r,t) 

x Idvl dw \v — w\p(r,v,t)p(r + s,w,t) . (10) 



w<v 



The finally resulting relation is, with 8v — v — V, 
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5p{r, v, t) « p (r, v, t)r (3 Sv/Q - 5v 3 /B 2 )/2 d r 
+ T (l-p)x(r + l,t)B(v) 

- P Q (r,v,t)T 6v/Q(l-p)x(r + l,t)fdv5vB(v) 

- P (r, v, t)r /(2Q 2 )(5v 2 - 6)(1 - p) X (r + I, t) 
xJdv5v 2 B{v) . (11) 

With p(r,v,t) ~ po(r,v,t) + 8p(r,v,t), it follows that the quantities p(r,t), V(r,t), and 
0(r, t), which are taken into account by the Gaussian approximation Po(v;r, t), are not 
corrected by Sp(r, v, t). However, for the third central velocity moment V we obtain, instead 
of T w 0, 

r = -3v^e/[(i - p)px] d r & + s& a r v + s 2 e/2 a 2 v 
-3sV^e/2d r e-3s 2 V^e/4d 2 e, (12) 

which becomes different from zero in inhomogeneous traffic situations. This causes the 
additional contribution 

+ ps/(6v / ^e)9 r r + ps 2 /(i2v / ^e)9 r 2 r (13) 

to the velocity equation (§) and the extra term 

- T[p + s d r p + s 2 /2 d 2 rP - ps/V^&d r V 

-ps 2 /(2v^e) d 2 V] - ps/2 d r T - ps 2 /A d 2 r T (14) 

to the variance equation @. Together with the continuity equation @, the resulting equa- 
tions are the desired macroscopic traffic equations for high densities. The related instability 
diagram is depicted in Fig. 4 and indicates two different kinds of instabilities. 

In summary, we have found several significant results, which are not sensitive to the 
particular choice of the parameters or to different variations of the model: (i) A conistent 
traffic model needs to take into account vehicular space requirements as well as Navier-Stokes 
terms in order to allow a realistic description of traffic instabilities, (ii) Treating vehicles in 
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a point-like manner, the stability of equilibrium traffic flow with respect to perturbations 
of small wave lengths is not correctly described, even if Navier-Stokes terms are included, 
(hi) This problem vanishes when corrections due to the finite space requirements of vehicles 
are considered. That is, the macroscopic traffic equations for low densities do not allow an 
extrapolation to the traffic dynamics at high densities, (iv) The Navier-Stokes terms are 
responsible for a subdivision of the instability region into separate areas. These belong to 
different eigenvalues, (v) Whereas the instability diagrams of traffic models, which consist 
of a density and a velocity equation only, typically show two relevant humps @H, the two 
additional (narrow) humps of the above model are related to the dynamic variance equation. 
Therefore, the dynamic variance equation gives rise to a new kind of traffic instability. 

A more detained discussion of the applied relaxation time approximation as well as of 
equations, simulations, and results for the non-linear regime of traffic dynamics will be given 



in a forth-coming paper [14 
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Fig. 1: The instability diagram of the 
Navier-Stokes-like traffic model for point-like vehicles (s = I = l® = 0) indicates that traffic 
flow would be unstable above a certain critical density, surprisingly even at large wave numbers k. 




Fig. 2: Same as Fig. 1, but for macroscopic 
traffic equations calculated up to Euler order only, thereby neglecting the deformation of the 
velocity distribution in inhomogeneous traffic situations. 
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Fig. 3: Same as Fig. 2, but with consider- 
lliaN: , / = 0.8s • V). In agreement with empirical 



ation of vehicular space requirements (Iq = 1/p 
findings the results predict that equilibrium traffic flow is stable up to 12 vehicles per kilometer 
and lane, at extreme densities, and at high wave numbers (small wave lengths). 
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Fig. 4: Same as Fig. 3, but for macroscopic 
traffic equations including vehicular space requirement and Navier-Stokes corrections in order to 
obtain valid results. The instability diagram is now divided into two separate humps in each 
half-plane, indicating two different kinds of instability. 
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